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^ ! 1 Introduction 

Random matrix theory has led to the discovery of novel matrix models and 
novel statistical distributions, which are defined by means of Fredholm determi- 
nants and which, in many cases, satisfy nonlinear ordinary or partial differential 
■ equations. A crucial observation is that these matrix integrals, upon appropri- 

ate deformation by means of exponentials containing one or several series of 
time parameters, satisfy (i) integrable equations and (ii) Virasoro constraints 
with respect to these time parameters. Most of the time, such matrix integrals 
can be written — by expressing the integrand in "polar coordinates" — as a 
multiple integral, which then can be expressed in terms of the determinant of 
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a moment matrix; this may be a moment matrix with regard to one or several 
weights. The extra time parameters are added in such a way that each weight 
has its own exponential time deformation. 

The main point is to show that this determinant satisfies (i) and (ii) . These 
features turn out to be extremely robust! The purpose of the present paper is to 
show point (i) in great generality, which is the determinant of moment matrices 
associated with one or several weights and defined on various different domains, 
satisfies the multi- component KP hierarchy with regard to the time parameters. 
This is a very general class of integrable equations. 

This determinant will turn out to be the r-function of this integrable hierar- 
chy; this r-function with appropriate shifts of the deformation variables will be 
expressed in terms of the "orthogonal polynomials" defined by the weights and 
their Cauchy transform. We list below a number of examples having their origin 
in Hermitian random matrix theory, in random matrices coupled in a chain, in 
random permutations and in Dyson Brownian motions (non-intersecting Brow- 
nian motions) on K leaving from the origin, where some paths are forced to end 
up at one point and others at another point, etc. . . These examples will then be 
discussed in detail in Section [7J 

• GUE: orthogonal polynomials. 

[ A!(z)TT e^i**** p(zi)dz t = detf [ z l+ i e^=± tkzk p{z)dz) 
n\J E n Vr ' 



0<i,j<n-l 



Coupled random matrices / Dyson Brownian motions: bi-orthogonal 
polynomials. 

■i [[ A n (x)A n (y)f[e^^- s ^p(x e ,y e )dx e dy e 

= det ( // x i y j e^=^ tkxk - SkV ^p(x,y)dxdy 



0<i,j<n-l 



• Longest increasing subsequences in random permutations: orthogonal 
polynomials on S 1 . 



1/ iA„(z)i 2 nf^ 



27r\/— lze 



V/s 1 2"7rV^lz / 0<ij<n-l 

mi + TO2 non-intersecting Brownian motions on M leaving from and 
mi paths forced to end up at ±o: multiple orthogonal polynomials on 



1 



A mi _)_ m2 (x, y) 
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'"1 2 \ 

A mi (a;) [| e-^ +axe e^=^ tk - Sk ^dx e 
v e=i / 

f m 2 2 _ \ 

v 1=1 J 



det 



\ \Je J °<'<™2-l / 

' < j < m x + m 2 - 1 / 



< « < mj - 1 

< 3 < "H + m 2 - 1 



• X)a=i m « = y^fl— 1 non-intersecting Brownian motions on R 7 with 
m a paths starting at a a el and paths forced to end up at bp: 
mixed multiple orthogonal polynomials (mixed mops) on R. 

A moment matrix for several weights: Define two sets of weights 

ip 1 (x),...,ip q (x) and tpx(y), . . . ,ip p (y), with x, y G R, 

and deformed weights depending on time parameters s a — (s a i,s a 2, . . .) (1 < 
a < g) and tp — {tpi, tp2, . . .) (1 < (3 < p), denoted by 

ip~ s (x) := ip a (x)e~^*'= lS °' kxk and <Pp(y) ■= 1 0fj(y)e^= 1 t0hV . 

That is, each weight goes with its own set of times. For each set of positive 
integer^ 

m = (mi, . . . , m q ) 1 n = (m, . . . , n p ) with \m\ = \n\, 

consider the determinant of a moment matrix T mn , composed of blocks and of 
size \m\ = \n\, with regard to a (not necessarily symmetric) inner product (■ | •) 

Trnn (^1 , ■ ■ ■ 3 Sq ) ^1 3 ■ ■ ■ 3 ^p) 



det T„ 



(T 7 



det 



11 

run 



rplp 



rpqp 
mn 



det 



\ mn 

f((x^ s (x) |yVi(y)}) < i<mi 



s (x) I y j <pl(y)))o<*< mq 

\ 0<j<»i 



l \ m \ = J2a=l m a and M = Z)fl=i n /3- 



0<j<n p 



0< 
IX J 



(1) 
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A typical inner product to keep in mind is 

(/(*) \g(y)) = II f(x)g(y)dn(x, y ), (2) 



where fi = [i(x, y) is a fixed measure on Mr, perhaps having support on a line 
or curve. 

From moment matrices to polynomials and their Cauchy transforms: 

I. Then, for 1 < /3, (3' < p, the following expressions are polynomials (with 
coefficients depending on s and iJH 

e ^\~n a ,-i Tm < n + e e- e P>( t P' ~ [fZ_Jj ._ n(P>P')(^ J of degree < 
e/9/3'W* — - y m „ (zj, jfor^^ 

(3) 

satisfying, for each /?, the following orthogonality conditions 



II. In the same way, the following expressions are polynomials (depending on s 
and t) 

e? a (n, m) z m ^ T "— + = P n *tf.a)( 2 ) of degree < mQ (5) 

7m n 

satisfying, for each /?, the orthogonality relations 

' (t | ) - » f » { ^ |1 ^/-C-'i 1 

/ 9 



X)^' tt) (^a'(») 



\a=l 



(6) 



III. The following expressions are Cauchy transforms of the polynomials ob- 
tained in II: 



introduce the notation [a] := (a, . . .) for a <E C. Only shifted times will be 

made explicit in the r-functions; i.e., T mn (t( — [z~ x ]) means that r mn still depends on all 
time parameters, but the variable t( only gets shifted. Moreover, here and below all the 
expressions £ a /3( Tl ); e ap{ n > m )i etc... all equal ±1 and will be given later. Throughout the 
paper, we use the standard notation e\ = (1, 0, 0, . . .), = (0, 1, 0, . . .). 



4 



^ Wt, + [*-*]) = ^ P ^ ){x)i) - S{x) 



z-y ! 

T " m \a=l Z ^ 

(7) 

IV. Similarly, the following expressions are Cauchy transforms of the polyno- 
mials obtained in I: 



, , -m a -l Tm+e a< n+e (s a - [z X ]) / 1p a s (x) 

e af 3{m,n)z a = ( 



v 



Q<£fHvH>(y) 



p'=i i 

(8) 

The statements I, II, III and IV summarize sections 1, 2 and 3. As will appear 
in Section 2, the polynomials appearing in (I) are called Type II i t mixed 

multiple orthogonal polynomials, whereas those appearing in (II) Type 1 1 t 
mixed multiple orthogonal polynomials. These were introduced by E. Daems 
and A. Kuijlaars [9], in the context of non-intersecting Brownian motions; they 
are a generalization of multiple orthogonal polynomials, where instead of one set 
of weights, there are two sets (the classical orthogonal polynomials correspond 
to one set with one element). They were introduced and studied by Aptekarev, 
Bleher, Geronimo, Kuijlaars, Van Assche [H El [T5l [8] . Around the same time, 
they were introduced by Adler-van Moerbeke in the context of band matrices 
and vertex operator solutions to the KP hierarchy [2] . In [7[ [5] , they were used 
in the context of non-intersecting Brownian motions and random matrices with 
external source. 

The (p + g)-KP hierarchy: Define two matrices W mn {z) and W^ni 2 ) °f 
size p + q, whose entries are given by ratios of determinants r mn of moment 
matrices as above, but with appropriately shifted t and s parameters. They 
turn out to be the wave and dual wave matrices for the (p + q) -KP hierarchy. 
It is remarkable that, upon setting all t and s parameters equal to zero, the 
matrix W mn (z) below is precisely the Riemann-Hilbert matrix characterizing 
the mixed multiple orthogonal polynomials! Similarly W^ ln (z) at t = s = 
satisfies the Riemann-Hilbert problem characterizing alternately the "dual" 
multiple orthogonal polynomials or the inverse transpose matrix of W mn (z) at 
t = s — 0. The Riemann-Hilbert matrix for the multiple-orthogonal polynomials 
has been defined in Daems-Kuijlaars [5], which is a far generalization of the 
Riemann-Hilbert matrix of Fokas-Its-Kitaev [TT] and Dcift-Zhou [TU]. Using 
identities as in I to IV, the two left blocks of W mn and the two right blocks of 
W^ n are mixed multiple orthogonal polynomials, and the remaining blocks are 
Cauchy transforms of such polynomials; for explicit expressions, see Section 
The matrix W m n(z) is defined by 
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W mn {z) diag (e~ ^ tlkzk , . . . , e~ £~ , e" s i^ fc , . . . , e ~ Ef 

3''"' ^"T^ ^ 13 00 ... (e a0 (m,n) 



l<0'<p l<o<g 

/ m -e n ,^.^' r [^'l) A / ^ m + ea -e Q ,,„(^-[^ 1 ]) 5 ,_!_ m \ 

■vmm.n) Z P £^/^ 1771 J 2 aa 

, y Tm " y !<«<<, ^ " aK ' y !<„'<„ 

\ l</3<p l<Q<q / 

(9) 



with inverse transpose matrix given by 

W* mn {z) diag (e E ~ tlkZ " , . . . , ***** , s ^ z \ . . . : 



l<P'<p \ mn / l</3<p 

l</3<p l<o< 9 

_ A ( ^ m+eQ _^,„( aQ /+^- 1 ]) 5 ,_ 1+m \ 

c Pa (n, m ) ^ * j ^ (^WM r= »«■-' «'J ^ 

\ l</3<p l<a'<9 

(10) 

The matrices W mn (z) and (z) satisfy the bilinear identities which char- 

acterize the t -function of the (p + q)-KP hierarchy 

W mn {z; s, t)W: n , n , 0; s*, t* ) T dz = 0, (11) 

for all m,n,m*,n* such that \m\ = |n|, \m*\ = |n*| and all s,t,s*,t* £ C°°. 
The integral above is taken along a small circle about z = oo; writing out the 
identity above componentwise and using the expressions ([9]) and (| 10[) for and 
W*, the bilinear identity (|11[) is equivalent to the single identity 



f (-l)^ ( " ) r,„,„_ C/3 (t /3 -^- 1 ])r m . i „, +C/3 (tJ + [ Z - 1 ])e E ? 0(t ' 3fc ~ t ^» ifc z ^-"J 2 dz= 
® (_l)"<*( m ) r m+ec , j „(s a ,-[z- 1 ])r m «_ e , ctira . ( S ; + [z" 1 ])e E ?° (i! °fc" s i* ! )ifc z m Q- m °- 2 <fe, 



where \m*\ = \n*\ + 1 and |m| = |n| — 1 and 

cr a (m) = (nw - m*,) and ^(n) = {n? - n* 0/ ). 

a'=l (3' = 1 

It remains an open problem to have a clear understanding of why the W mn (z; s, t)- 
matrix above, evaluated at t = s = 0, coincides with the Riemann-Hilbcrt 
matrix for the mixed multiple orthogonal polynomials. 
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PDE's for the determinant of moment matrices: Upon actually comput- 
ing the residues in the contour integrals above, the functions r mn , with |m| = \n\, 
satisfy the following PDE's expressed in terms of the Hirota symbol: 



d 2 

''dSaJ+idSa',! 
d 2 

d 2 



In Tmn 
111 Tmn 



hi: 



In- 



>SV+2c5 (3/3 , (dtp)T mt n+ep—e.pi ° Tm.n+e^i — ep 
Sl+28 aa , (9s a Pm+e a i —e a ,n ° T m+e a —e a i ,n 
S^Otp )'7"m+e a ,n+ep ° 7"m — e a ,n— ep 
£*({d Sa )7Vn — e Q ,n — ep T m+e a ,n+ep ■ (12) 



Whereas the formulae above have in their right hand side different r m „'s, one can 
combine these relations to yield PDE's in a single r mn ; so, these are PDE's for 
the determinant of the moment matrix {T]). In particular, one finds the following 
( p ~^ q ) PDE's, which play a fundamental role in chains of random matrices and 
in the transition probabilities for critical infinite-dimensional diffusions: 



d 



In T mn 



d 



Ql i nT 



■ lnr„ 



d 



ds a ,!ds all 111 'mn 
dtp, 2 d ScCtl 111 Tmn 



ds a ,i 



In: 






a 2 


In T~mn 


dtp 






d 2 


In T~mn 


dtp 






d 2 


In 7"mn 


ds a , 


>2 9Sc.,l 




d 2 


In T mn 


ds a , 


,l3So,l 




d 2 


In T m n 


ds c 






d 2 


In T mn 







= o, 



= o, 



= 0. 



2 Tau functions and mixed multiple orthogonal 
polynomials 

Following [9] we introduce the notion of mixed multiple orthogonal polyno- 
mials (mixed mops), with regard to two sets of weights {</?i,<£>2, ■ ■ • ,<^p} and 
{•01,^2, ■ ■ ■ ,tpg}- 



e For a given polynomial p{ti, <2, ■ • • )> the Hirota symbol between functions / = f(ti, t2, ■ . .) 
and g = <?(ti, *2, ■ ■ •) is defined by: 



P( 



a a 
at? at? 



.)fog :=p( 



a a 



■)/(* + y)9(t - y) 



y = 



We also need the elementary Schur polynomials Si, defined by e^i t ^ z := ^2 k>0 Sk(t)z k 
for I > and Sg(t) = for I < 0; moreover, set 

" ; eK dti'2dt 2 3dt 3 
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Definition 2.1 Let A\, A 2 , . . . , A p be p polynomials in y and set 

Q(y) := A^ip^y) + A 2 (y)ip 2 (y) + ■■■ + A p {y)tp p {y). 

Type I For a <G {1,2, ... ,q} the polynomials A\, A 2 , . . . , A p are said to be 
Type I normalized with respect to i\) a , denoted Type I | ; */deg(j4^) < np for 
= 1, . . . ,p and Q satisfies the following orthogonality conditions 

(x l ip a ,(x) | Q(y)) = <W<5 4 ,m Q -i, i = 0,. ..,m al - 1, 1 < a' < q. (13) 

Type II For (3 G {!,..., p} the polynomials A\, A 2 ,...,A p are said to be 
Type II normalized with respect to tpp, denoted Type II ^ , if Ap is monic of 

degree np and deg( Api) < np> for 1 < f3' < p, with 0' ^ 0, and Q satisfies the 
following orthogonality conditions 



(s'lM*) \ Q(y)) = 0, 



2 = 0,..., m a — 1, 1 < a < q. 



(14) 



In both cases, the polynomials A\, . . . ,A p are called multiple orthogonal poly- 
nomials of mixed type, or mixed mops for brevity. 



Proposition 2.2 For [3 = 1, . . . ,p, let 

Q%Uv) ■■= Q%P(v)<p\(v) + ■■■ + Q^ p) (yH(y), 

where Q ( £f\ with 1 < 0, [3' < p are the polynomials, defined by 
jni j mn {tp- [z- 1 ]) 



(15) 



Qm.n (•£■) 



Q ( if\z) := epp,{n)z n ? 



0'+0, 



and 



(_l)« / 3' + l+n / 3' + 2 + ---+™^ + l > 01 ^ 

(_l)"^+i+ n f <+2 + ---+n,5' if < 0' . 



(16) 



(17) 



Then Qmn\y), • • • , Qmn'(y) are Type II | t mixed mops. 



Proof For j = 0, 1, 2, . . . and = 1, 2, . . . ,p we define a column vector Ctj 



of size \m\ by 



._ 



' ((xHi s (x) y^p(y)))^. 

v \ / / o<ii <mi 

((a^'C*) |^(y) 



((^V'q s {x) y j <pp(y) 



0<i 2 <m 2 



0<i„<m„ / 



(18) 
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When its size is important (see the proof of Proposition 12.3ft we write Cj (m) 
for (fT8)) . Notice that the moment matrix T mn , defined in can be expressed 



in terms of the columns Cj , and so 



det 



\pQi Cfj ■••) Q)i •••) C« p -l) • (19) 



For future use, let us point out that the dependence of T mn on the t variables is 
as follows: 

r mn (h) = detail), CUh), C^iti), Cl Cl C^_ x ), 
r mn (i 2 ) - det(C \ C{, Ci^, C*(t 2 ), C?(i a ), C£ A 



Tmn(tn p ) — det(C , C\ , C^ 1 _ 1 , C Q , Cf , C£ p _ j (t„ p ) ) . 

(20) 

Since 

p p / n^z-l \ 

= E Qmn'\y)<Pp'iy) = E + E 4?** Mv)> 

/3' = 1 /3' = 1 \ j=0 I 



(21) 

the orthogonality conditions (JT4J) for Q = Q™ can be written as the linear 
system 



E E A 0P' C 3 =~ C n^ 
f}'=l j=Q 

of |m| equations, in the |n| (= |m|) unknowns A^,, where 1 < /?' < p and < 
j < np>-\. If we order these unknowns as follows: A^, Ag X , . . . , , A°p 2 , 

Ap 2 , ■ ■ ■ , A n ^ v , then this linear system has precisely r m „ as determinant, in 
view of (fT9]) . Since r mn 7^ 0, generically, we have by Cramer's rule, 

A^, = — ^ f— ? (22) 

' mn 

Substituted in (I2ip this yields an explicit expression for the Type II 1 mixed 

mops Qmn\y), . . . , Qmn\y)- 

In order to connect these polynomials with the tau functions r m „ we first 
expand T mn (tp — [z^ 1 ]) using ([20]) . Thus, we need to compute Cj(tp — [z^ 1 ]), 
which we claim to be given by 

Cf(tp - [z- 1 ]) = C?(tp) - z- l C? +1 {t p ) =Cf- z- l C? +v (23) 

where the last equality is the notational simplification agreed upon. To prove 
the first equality in (|23|) . which is an equality of formal series in z~ 1 , let us write 
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a typical entry of the column vector Cjitp) with its explicit time-dependence 
on tp, 



where 1 < a < q and < i < m a . The following trivial identity will be used 
over and over again in this paper 

e -Er'T = i_ l: . (24) 

In view of the latter, the same entry of Cj(tp — ) ( as above) is given by 

(x^- s {x) \y%(y)(l-l)) = (^( X ) \y^ P {y))-- z {x^- s {x) \yi +1 ^{y)) 

which proves (|23|) . Using the fact that the determinant is a skew-symmetric 
multilinear function of its columns, which vanishes when two columns are equal, 
it follows from (gDJ), (HD and ^ that 

— act^OQ, . . . , o„ (3 i _ 1 , zuq , . . . , o n(3 , u , . . . , o nj) _ 1 ; 



J=0 

™/3 



E 7 j Apttr 1 r 1 r 13 ~r< 3 r" 3 r 13 r^ +1 r p "i 



i=0 

np 
.7=0 

In (*) it is understood that all the columns between —Cj +1 and —C^ come with 
negative signs and no others; this notation shall be used freely in the sequel, 
without further mention. 

For Qmif ' with (3 ^ j3' we also need to keep track of signs and of shifts 
in the first index of the tau function, as is seen in the following computation, 
where we suppose that f3 < (3': 

z n e'-iT m , n+ep - ef3 ,(t , - [z- 1 ]) 

- deter 1 r 13 ri 3 ?c 0, - rZ zr 13 ' -r 13 ' r p ' +1 r p ) 
= ^2 zJ det(Co, . . . , Cf a , . . . , Cf 3 _ 1 , —C? +1 , . . . , -C^,_ i , +1 , . . . , C T ^ p _ 1 ) 
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e m {n) Y, * J det(C 1 , . . . , C^_ 1 ,Cl + \ ...,cf_ lt cf +1 ,..., C^_ x ) 



3=0 



3=0 

= ei )f3l {n)T mn Q { if ) {z). 

The sign eppi{n) which we introduced when moving the column C@ to the right 
is given by (— in agreement with (I17p . When (3 > f3' the column 
is moved to the left, which yields a sign eppi{n) = — (— V) n 0'+^ h ™' 3 , as is 
easily checked. □ 

The tau functions r mn also lead to Type I normalized mixed mops, as given 
in the following proposition. 

Proposition 2.3 For a = 1, . . . , q, let 

PiuHy) ■= Ph'HvViiv) + ■■■ + PL a /HyV P (y), (25) 

where Pmn^ are the polynomials, defined by 

P^fXz) := e a0 (m,n)z^ Tm - e - n - e ^- [z ~ 1] \ (26) 
with leading sign 

e al3 (m,n) = (27) 
XTien Pi^ 1 ^?/), . . . , Pmn P \y) are Type I 1 _ mixed mops. 
Proof Letting 

p^liv) = E ^(i/HM = EE B ^ y 3 ( 28 ) 

/3=1 /3=1 j=0 

the orthogonality conditions (H2J) for Q = ^mn can be written as the linear 
system 

p mp—l 

E E B U c i= E ^ 

0=1 j=0 

where denotes the column vector of size |m| with a 1 at position m a of the 
a-th block (so at position mi + m-i + ■ • ■ + m a ), and zeros elsewhere. Cramer's 
rule now yields 

j det(Co, C{, . . . , Cj-i, Em a i Cj + x, ■ ■ ■ , C^_ 1 ) 



"a/3 



1)j+ i- n , ^t{DlDl...,D^ 11 D^Dl +1 ,... 1 Dl p _ 1 



e Q/3 (m, n)(-l) 
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where the last line was obtained by expanding the determinants along the 
column, e a p(m,n) is given by ((27)) and Dlj, is the column vector UJ. with its 
(mi + • ■ • + m Q )-th entry removed, i.e., := C2(m — e a ). This yields ex- 
plicit expressions for the Type 1 1 _ mixed mops. To connect them with tau 

functions, we notice on the one hand that the columns D7 appear in the ma- 
trices which define the tau functions T m _ ea and on the other hand that these 
columns behave in the same way (|2"3")l as under shifts. Therefore we can 
compute, as before 

2 j Tm — e a .n — ep(tp ' % ]) 

np—l 

- V zidetlD 1 D 13 -O p -D* 3 D p+1 D p ) 

3=0 

= E (-l)"^^ 1 ^ defcpj, . . . , Df_ l3 nj, Df +1 , Dl^) 
3=0 

riff-l 

= e a p(m,n) ^ z J B 3 a/3 T mn 

3=0 

= e a p (m, n) F I ( „ Q I j /5) (z) r TOn . 

□ 



3 Cauchy transforms 

We now show that certain shifts of the tau function, appearing in the Riemann- 
Hilbert matrix of [9], are (formal) Cauchy transforms. For a function F and a 
weight "0, define its Cauchy transform as 

C*G{z)~(^- \ G(y))=f2-^ I (x i ^(x) \G(y)), (29) 

i.e., our Cauchy transforms will be formal in the sense that we always think of z 
as being large, and this is precisely how it will be used. The first type of Cauchy 
transforms which we are interested in are given in the following proposition. 

Proposition 3.1 For a = 1, . . . , q and (3 — 1, . . . ,p, the Cauchy transforms 

of Qml(y) = Qmn\y)tp\(y) H 1- Qmn\y)Vp{y)i with respect to ip~ s can be 

expressed in terms of tau functions as follows. 

C+-.Q&{z) = e a ^n)z—i (30) 
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Proof The proof is based on an investigation of the moment matrix by row. 
Therefore we define, for a = 1, . . . , q and for z = 0,1,2,... the row R l a of size \n\ 

by 



0<i P <n P 



When its size is important we write R l a {n) for R l a (n). The moment matrix T„ 
can now be expressed in terms of the rows R l al and so 



= dot 



(31) 



The tau function which we need to compute is T m+ea ^ n+ef} (s a — so 
throughout the proof, R l a , stands for R r a i(n + ep) for all 1 < a' < q and i — 
0,1,2,.... Notice that the only rows which depend on the time variables s a = 
(s a i, s a 2, ■ ■ ■) are the rows R l a . Recall the dependence of ip~ s on s a as follows 

i/j~ s (x) = ip a {x)e- E "=i SakX " , 

so that, according to the identity (HH)) , when s a gets replaced by s a — [z^ 1 ], 
then ip~ s (x) gets replaced by ^a S { x ) f 1 + — + + ■ • • ) . It follows that 



KiSa-l*- 1 ]) 



X^- S { X ) ( 1 + - + ^ + 



where we introduce the convenient abbreviation n'p, := npi + 8/3/31 — (n + ep)p< ■ 
Notice that 

for < i < m a — 1; we stop at m a — 1 because the highest index i for which R l a 
appears in T m + eQ in + ej3 is i = m a . By recursively applying this formula we get 
that for < i < m a — 1 

R l a (s a — [z^ 1 ]) = R l a (s a ) + lin. comb, of lower rows. 

This leads to the first equality in 



Z ^~m-\-e a ,n+ep {^a \_Z ]) 



(32) 
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z -m a -l det 



R. 



m a — 1 



V 



R, 



m„ — l 



det 



R a {z) 



V C 1 / 



For the second equality, in which we have put 

R a (z) :-- [ (\^— kVi(y) 



0<ji<ni 



\y jp <Pp(y) 



it suffices to show that 

R™ a (s a - [z^ 1 ]) = z m « +1 R a (z) + lin. comb, of higher rows R l a . 
To do this, compare a typical entry of R™ a (s a — [z^ 1 ]), to wit 

(z m "V>r to (i + f + + ■ • •) I y'My) ) 
= ^(^ s w(© mQ + (i)" l ° +1 + ---) IvVMv) 

with the corresponding typical entry of z ma+1 R a (z) , to wit 



0<j p <n' p 

(33) 
(34) 



(fc s W (i + f + © 2 + -") |yV^(y))- 



It leads to the following explicit expression for (|34|) : 

m a — 1 



and hence to the proof of the second equality in 

In order to make the connection with mixed mops, we introduce for f3' = 
1, . . . f p the row (z) of size \n + ep\ = \n\ + 1 which has zeroes everywhere, 
except in its /3'-th block, nameljfl 



sf(v) = (0...0 (^') c 



0...0 



7 Recall that n'„, = npi + Sppi 
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Notice that with this notation, definition of R a (z) can be rewritten as 



It suggests the introduction of the following polynomials (in y) 

( Ri \ 



(35) 



S^(y) := dot 



Dm a — 1 

7?° 



\ rT^ 1 J 



(36) 



Expanding this determinant along its (mi + ■ • • + m a + l)-th row, which is the 
(unique) row that contains y, it is clear that if /?' ^ (3, then degSf^(y) < n^, = 
np<. In view of (f3Tj) we also hav^l 



Moreover, for any a' — 1, . . . , q and i = 0, . . . , to q ' — 1, we have by linearity of 
the determinant 



I3' = l 

( 



det 



E>m„ — 1 

/3'=1 



(37) 



3 See 1271 1 for the definition of e a fj(m,n). 
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det 



r>7Tl a — 1 



0<j<n', 



det 



/ 



jym a — l 

K> 



V K^ 1 / 



which is zero, since the latter matrix has two identical rows (i < m a i). This 
shows that 

e a p(m,n) QP 



e a/3 (m, n) gl 

~ ^apKV) 



are type II mixed mops, normalized with respect to (pg. It follows from Propo- 
sition 12.21 that 

-'~ ^-mnQtffHy), (38) 

m) (y)<pt(y) + --- + Q ( £h p) (yW(y)> it 



for anjO a=l,...,q. Since Qml(y) 



follows from ([32j) , j35]), (J38|) and (20|) that 



,/3'=l 



= e a/J (m,n) C - 3 £ Qj&f' 5 

= e a/J (m,n)C^-.Q<&(*) (39) 

This finishes the proof. □ 

Observe our proof shows, as a byproduct, that each Qmn(y) is expressible nat- 
urally as a determinant, like in the classical case, namely 

/ *? \ 



Qiil ( = e^^Kn) 



det 



/3'=1 



7?° 



V R 



m„ — 1 



(40) 



9 The formulas for the different values of a are all the same, up to a sign, as they amount 
to changing the location of a row in the evaluation of a determinant. 



1G 



We now get to the second type of Cauchy transforms which correspond to the 
Type I mops Pm2(y)- 

Proposition 3.2 For 1 < a, a' < q the Cauchy transforms of Pmn {y) = 

Pmn'^iy) fiiy) + • • • + Pmn' P \y) ( p i p (y) with respect to ip~ s can be expressed in 
terms of tau functions as follows: 



j/>~ ran V') 



7~mn (-5 a Z ] ) 



' ra-\-e a —e 



,,,fi(s a - [z x ]) 



(41) 
a' 4&2) 



Proof Up to a relabeling of the indices, the shifted tau functions in question 
were already expressed as polynomials in the previous proof. Let us show how 
this leads to a quick proof of (|4ip . Shifting the m a and np indices down by 1, 
it follows from (|32|) that 

/ R°i \ 



"T mn {s a - [z V \) = det 



H a a 

R a (z) 

7?° 



while the orthogonality relations ([57)1 become 



I3' = l 



= det 



V K^ 1 J 



ota' u z,m Q 



(43) 



Since degS^a < rig/ for f3' = 1, . . . ,p this means that the polynomials 



- s ia(y) ■ 



are type I mixed mops, normalized with respect to ip a s , so they coincide ac- 
cording to Proposition 12.31 with the polynomials Pmrl^iy), ■ ■ ■ ,Pmrf\y)- We 
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conclude, as in that 

— m a T mn( s a ~ [ z ]) 



1 



/9' = 1 



Similarly, one obtains (|42|) from (|30p by shifting m„ and down by 1; the sign 
in this case is determined (for a' < a) from the right hand side of (|43|) now 
taking the form 



det 



n 'a' + l 
nm a -l 

K» 

Ft 



V R?^ 1 J 



s a ' a (m) det 



( R i \ 



J^ 7n a — 1 

7?° 



V iC* -1 J 



□ 



4 Duality 

By interchanging the roles of the weights ip l a with the weights ipp S we obtain 
Type 1 1 mixed mops and Type II i _ mixed mops, expressed in terms of tau 

functions, leading to a duality. As a general rule, in order to dualize a formula 
one does the following exchanges 

q «-> p, m <-> n, ijj ^ ip, s <-> — t, x <-> y. (44) 

At the level of the indices, duality amounts to 

a •*-> /3, i j. (45) 

As for the mixed mops which we have constructed, they will correspond to new 
mixed mops for which we will use the same letter, but adding a star. Thus, 

p*(f3) p(a,/9) ^ p*(/9,a) q(/3) ^ Q*(a) n(W) 
mn nm s mn nm ' ^mn ^nm J ^mn ^nm 

(46) 

What happens to the tau functions T mrl ? To see this, pick a typical shifted tau 
function r m + ea ,-e ,,n(sa — an( i make its dependence on the weights and 
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on all times explicit, writing T m + eot - e , lH (s — [z^ 1 ] e a ,t;ip,ip). According to 
the above rule it becomes r ra + e/3 _ e(3 , , m (— t— [z^ 1 ] ep, —s; ip, tp) which is equal to 
Tm,n+ef3-egi ( s i t + [ z _1 ] G/3', i/>, <p), since transposing the moment matrix has no 
effect on the determinant, while it permutes the indices in the tau function, it 
permutes the time-dependence (with signs) and it permutes the weights. Thus, 



Tmn(tp-[z *]) *->■ T mn (s a +[z x ]) 

7~m— e a ,n — (^/3 \_Z J) <— > T~ m _ £a n _ e ^ (s Q ~t~ ]) 
T m+e„— e ,,n( s Q ~ [ 2 1 ] ) 7 "m,n+e |3 — e^; (^/3 + [ z 

and so on. Dualizing Propositions l2.2l and l2.3i we get the following proposition. 



Proposition 4.1 For a = 1, . . . , q and (3 = 1, . . . ,p, let 

P*%\x) := P^ix^ix) + ■■■ + P^< q) (x)ip-°(x), 
Qlt\x) := Ql a ,' 1 \x)^ s (x) + --- + Qly\x)^(x), 

where Pnm' a ^ and Qnm' a ^ are the polynomials, defined by 



(47) 



P*J£^(z) := ep a (n,m)z m «- l Tm - e *' n - ei3iSa+ ^ ^\ (48) 



and 



n *(a,a), s n m a T mn(s a + [z *] ) 

tynm \Z j • — Z 

T~mn 

Q±' a \z) := e aa ,(m)z m «'- l Tm+e «- e <"' niSa ' + [ *~ J) , a'^ a. 

(49) 

T/ien Pnm (x), ■ ■■ i Prim (x) are Type I i mixed mops, luMe Qnm (a), ■ • ■ , Qnm' (#) 
are Type II i _ mixed mops. 

□ 

Dualizing Definition (|2"9"|) we get the following definition for the dual Cauchy 
transform: for any function F and a weight ip we put 



C\/m:) := // ^F(x)d^x,y) = (F>,,) 



<p(y) 



(50) 



z-y 

If we dualize now Propositions l3 . II and !3 . 2\ then we get the following proposition. 
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Proposition 4.2 For a = 1, . . . ,q and /3, /?' = 1, . . . ,p, the Cauchy transforms 

of Pnm \x) with respect to (fg, and of Qnm\x) with respect to ifg can be ex- 
pressed in terms of tau functions as follows: 

r* P*(£)M - T mn{tp + [z *]) 



' mn 

C;.P (1 ;.;V.< = :/)|;<) .^l-n, W + e,-e 1 ,fa+[.- 1 ]) i ^ ft 



and 

13 ''"mn 

□ 



5 The Riemann-Hilbert matrix and the bilinear 
identity 



Orthogonal polynomials were shown to be characterized by a Riemann-Hilbert 
problem in [TT] and [TUj- This was generalized by Daems and Kuijlaars to the 
case of mixed mops. According to [90 the corresponding Riemann-Hilbert 
matrix is given by the (p + q) x (p + q) matrix 



V 



\Wmn I l</3<p I ^tj,- s Wmn\ 1 



P 



!</3'<p 



l<a<(j 
!</5<p 



A> p("') 
^„/. — s r van 



i<p< P 

1</3'<P 



<P<P 
1<<V<(, 



l<a'<q 
l<a<q 



^p(m,n)- 



- 1 ) 



V 

whose inverse transpose matrix is given by 



l<a<q 
l<f3<p 



£ a 'c( m )- 



(r* p*(P')\ ( 

1</3<P 



l<a<q 
l<P<p 



- P 

r 7 



l<f}<p 
l<a<q 



( n *(a,a')\ 
\Wnm J i< a < q 



l+ea ,n+e»C*p + 



l</3'<p 
1</3<P 



l<o<g 
l<fi<p 



l<a <q 
e 0a (n,m)- 



e aC( t (m) - 



'Up to a factor diag(/ p , — 2ny/— ll q ) which we suppress. 
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We will obtain bilinear identities for these tau functions from an identity which 
is satisfied by the Riemann-Hilbcrt matrix and its adjoint. We define the wave 
matrix W mn (z) by Y mn (z)A(z), where A(z) is the diagonal matrbj"1 

A(z) :=diag(e* (tl '*\ ...,e* ( *» ,a) , e« (si ' z) , e« (s " z) ), 

with adjoint wave matrix Y^ ln (z)A~ 1 (z). In order to make the dependence 
on the time variables (s,t) explicit, we will write W mn (z; s,t) for W(z) and 
W^ n {z;s,t)ior W^ n (z). 

Theorem 5.1 The tau functions r mn satisfy the following bilinear identities 
that characterize the tau functions of the (p + q)-KP hierarchy (see 



W mn (z; s, t)W^ n , (z; s* ,t*) T dz = 0, 

which is equivalent to the single identity 
v , 

J2 f {-l) a0M r m ,n-e,(tf} - [z- 1 ])T m . in . +e0 (t}+ [z- 1 ])e« t "- t /"*^ n "- n /»- a dj = 

J2 I (-l) CT " (m) r m+e „,„(sa- [z- 1 ])r m ,_ ea ^(s* a + [z- 1 ]^*"-**^ z m «- m »- 2 dz 
where 

a (i 

<r a (m) = }^ (m a > - to*,) and 073(71) = } j (np> -n* gi ). (52) 

a' = l 0'=1 

and | m* | = | n* | + 1 and \m\ = \n\ — 1. 
Proof For the entry (/?',/?") of the product 

YA{Y*A~ 1 ) T = Ydiag^* 1 . . . , e^~ l >\ S ai —i>*\ . . . ,e^ s "-<'^)Y* T , 
we need to prove that 



f3=l J °° c t =l J °° 

(53) 

where it is understood that all polynomials P* go with starred times s* and 
t*. Also, the integral stands for (minus) the residue at infinity, and can be 
computed using the following formal residue identities, with f(z) — X^o ^*' 

' ' f(z)C^g{z)dz = (f(x)iP(x)\g(y)), (54) 



2ttV-1 

' c;f(z)g(z)dz = (/(*) W(y)g{y)) , (55) 



2W-1 



1 Throughout this section, we set £(t, z) := *fc 2 * 
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whose proof we defer until the end. Using this, and Definition (jlT]) of the 
functions Pnm\x), the left hand side in (f53|) becomes (up to a factor 2 7r -*/—!) 



13=1 (3=1 

= (i#£?(x) \Q%>(y)). 
Similarly, the right hand side in becomes (up to a factor 2tt\/— T) 



£(i#& a) (*)l&« - 0» 



o £0 Q -s* ,x) 



a=l 

= (p:C\ X )\Q^(y)). 



The three other identities are obtained in the same way. 

In terms of tau functions, it means that we have shown that for any to, n, m* ,n*,(3' 
and f3", with \m\ = \n\ and \m*\ — \n*\ the following bilinear identities hold: 

J2 i * T ™,n+e fjl -e l3 (t l 3- [z~ l ] )r m . ,n'+e^e „ [t% + [z~ X ] )^~**^ dz = 
(3=l Jo ° 

<f * T m+ecttn+e0l {S a - [z- 1 ] )T m *-e a ,n*- efjl , (s* + [z' 1 ] )e« (s ^ s ^ 2) dz, 
i >^ oc 



where 



* = e p ,p(n)ep„ (n*)z n ?- n v- 2+s ^' +s ^", 
4 = e al3 ,(m,n)e f3 , la (n*,m*)z m *o- m <*- 2 . 

For different values of j3' and /3" this yields the same identity, up to a relabeling 
of n and n*. Namely, replace in the bilinear identity n + epi by n and n* — epn by 
n* and multiply by (-l) n i+-+ 7 V (— l)""* -1 l_ri 3" to find the following symmetric 
expression for the identity, that is independent of (3' and (3": 

£ <f (-l)' j( "»V-e,(t r [2- 1 ])r m ., n . + e /J (iJ+[«- 1 ])e« t "- t 5'*^ n '»- n 5- a d? = 

where <t q (to) and crp{n) are given by (|52| . Notice that, due to the shift, one 
must have in this symmetric form that |m| = \n\ — 1 and \m*\ = |n*| + 1. The 
other three identities also yield the above identity, up to relabeling. 



22 



Finally, to prove (|54[) . compute 



\ i / °° 00 1 

1 9(y) ) dz= j> 22 a t z l 22 — (x^(x) I ) 

OO / OO \ 

= J2 a * 1 9(y) ) = ( E ^ ^ 1 9(3f) ) = (f(x)1>(x) I g(y) ) , 

i=0 \ i=0 I 

and similarly for (|55|) , completing the proof. □ 

6 Consequences of the bilinear identities 

In this section we will derive from the bilinear identities (|5ip a series of PDE's 
for the tau functions T rnn . In order to keep the formulas transparant we will use 
the following simplification in the notation. Recall that we have time variables 
s a = (s a i, s a2l ■ ■ ■) and tp = (tpi, tp 2 , ■ ■ •), where a = 1, . . . , q and /? = 1, . . . ,p. 
In the bilinear identities (p)Tj) we consider in each term a shift in t a , for a single 
a, or in sp, for a single /?; we will denote this t a or sp by v (so v is an infinite 
vector v ~ (v\,V2, ■■ •) and we assemble all the other r := p + q — I series 
of time variables in w — (wi, W2, ■ ■ ■ , w r ), where wi — (ton, W12, . ■ .) and so 
on. Moreover, precisely like in the bilinear identities we will want to consider an 
independent collection of all these variables, in fact we will consider here (v' , w') 
and (v" , w") besides (v, w). We use the Hirota symbol, which takes in our case 
the following form 

P{d v ,d w )FoG = P(d v <,d wl )F{v + v',w + w')G(v-v',w-w')\ . (56) 

The elementary Schur polynomials Se(v) are defined by 

00 

e^"*** =J2$k(v)z k , (57) 

fc=0 

for I > and Si(v) :— otherwise. In particular, if we put degree Vi := i, then 
50 = 1, Si(v)=vi, St(v) = V£ + degree £ in v\, ... ,vi- X . (58) 
We also use the standard notation 



d,, 



d 1 d 1 d 

dv\ ' 2 dv2 ' 3 9^3 



We first give an identity which will allow us to compute the formal residues 
which appear in (|51| in terms of derivatives of the tau function. 
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Lemma 6.1 For any n € Z we have the following formal residue identity 



F(v" + [z- 1 ] , w") G(v' - [z- 1 ] , w') e KSo(^-<)^ z» "' : 



= E Si-i-»(-2o) ^(<%) e ^»'< + ^ F(v, w) o G(v, w), 

j>0 

(59) 

where 

v' = v — a, v" = v + a, w' i =Wi — hi, w'/ — Wi + hi, 

a = (01,02,03, • ■ .)> h = (hi, b i2 , b i3 , . . .), 

for 1 < i < r. 

Proof The proof is an immediate, but tricky, consequence of Definition ([57| of 
the Schur functions and of the following two properties of the Hirota symbol: 



F{v+ [z- 1 ] ,w)G(v- [z- 1 ] ,w) =J2z- j S:j(d v )FoG, 

3=0 

F(v + a,w + b)G{v-a,w-b) = e ^=o(°^+E; =1 ^0%) p Q q 



□ 



Proposition 6.2 The bilinear equations imply, upon specialization, that the tau 
functions T mn , with \m\ — \n\ satisfy the following PDE's expressed in terms of 
the Hirota symbol: 

2 d 2 

r mn7T TjT hlT mn = Si+25 00 , [dtp )T m ,n+ep-e , ° T m ,n+e , -e p (60) 



1 dtp^+idtpi ,i 
d 2 

T mri^ n hlT mn = % + 2A . (9 s „)r m+e ,_ e n o r m+e e . n (61) 

C/Sa^+lOSa/,! 

"ia m In^mn = S e (dt )T m+ea . n+e0 o r m _ eQ! „_ e/3 (62) 



3 2 



i 



lnT mn . — Sf,(d Sa )T m —e a ,Ti—e0 ° Tm+e ,»+e|! (63) 



Equations A60\) resp. H61)) for f3' — (resp. for a' = a) yield a solution to the 
KP hierarchy in tp (resp. in s a ), while for (3' ^ ft and a' ^ a, h6U\) — i6S\) 
yields 



d T m,n+eg-e g / T m,n+e g / —en 

— — lnr„ m = '— (64) 

OT/3,lOT/3',l 



rii n 
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7"m+e / — e a ,n^"m+e a — e_,/ ,n 

mr m „ = 5 (65) 



ds a Ads a > i 



9 2 7"m+e Q ,ri+e« T Vre-e Q ,ra — e« /v./.\ 

lnr„ m = = ( 66 J 



ds a ,idtp 



ran 

d 2 



■In 



Wf — ln = a 3 i ( 67 ) 

01/3,1 T- m: „ +ef) ,_ e/3 ____l n7 - mn 
9 Tm-e a +e„/,n da a , 2 ds a , A n T > 

« — ln — = — & — , ( 68 ) 

GSot.l ^~7n — e n /-\-e a ,n q- or 111 T 



niit 



— Ill Tm+eQ,n+e ^ _ dtp^dsa,! ^ nTrnn {[)*)) 

dtp,! r m _ ecv;n _ e/3 —d-— lnr mn 



— In ^ = . (70) 

- lnr„ 



/a a 



It leads to the following ( p ~^ q ) PDE's for lnr mn involving not just one s a or tp, 
but a few of them 



tor,! \ m^k^ lnTmn I dtp 
, §1 mT 

O I ds ai2 ds a , 1 111 ' mn 



Id 2 \ 





d 2 


In T~mn 








d 2 


hi T~mn 








a 2 


hi T~mn 








a 2 


In T mn 


v 9s a / 


tl ds a ,l 




d 2 


In T mn 


[ ds c 


i,23t/3,l 




a 2 


In T mn 


\ ds c 





- ~ i — ^ I = 0, (71) 

= 0, (72) 



0. (73) 



Proof Let us denote for a = (ai, . . . , a 9 ) and b = (b\, . . . , b q ) by fi(a, 6) the 
differential operator 



= 1 \a' = l /3' = 1 ' 



Using Lemma (|6.ip . rewrite the bilinear identitjo ([ST]) : 

p 

- e/3 



E(-l) CT " (n) £5 n; _^ +1+fc (-26 /3 )S fe (^)e n ^ i, )r m ., n . +e)3 o r„ 



£=1 fc=0 
g oo 



£(_l)ov.(m) ^^ mQ „ m - +1+fe (-2a Q )^(9 s Je 0(a ' fc) r m ._ ecv ,„. or m+ea ,„ = 0. 



fc=0 



(75) 



12 Recall that in this form of the bilinear identity |m*| = |n*| + 1 and |m| = |n| — 1. 
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Note that all infinite vectors a a and bp can be chosen completely arbitrary. We 
set all components of a and b equal to zero, except 6/3^+1 = B 7^ (for some 
fixed P and £), and we set to* — m and n* — n = —2epi (for some fixed /?'). 
Then only the first term in (JT5J) survives, the signs 0/3(71) are all 1 (see (|52|l ) 
and, in view of (|58p. the identity ([75)1 becomes 



P OO 

= £ ^5l+fe-2V^(- 2 V')^(9 V ') eiJ77 '^+ rT -m.„ +e(3 ,, 
/3"=1 fc=0 

( 52 \ 2 

= -B ^ — 2SI+2&PPI \ptpPm,n+ep— 2e ff , Tm,n-ep + 7^ ^ Tm,n—e^i T m ,n-e ff , J + 0[B ). 

Expressing that the coefficient of B in this expression must vanish we get (|60p , 
upon relabeling n—epi i— > n and upon using the following property of the Hirota 
symbol, valid for / depending on (time-) variables s and t: 

■ F o F = 2F 2 — — lnF. (76) 



<9i9s dtds 

(j6T|) follows from © by duality, using P(-<9 S ) FoG = P(d s ) GoF. In order 
to obtain ([62]) we consider again (f75|) , with 6/3,^+1 = B 7^ and all other 
components of a and 6 equal to zero, but we set now n* = n and m—m* = — 2e a . 
Then (l75l) becomes 



= ^Sfc+i (-2bp)Sk (dp) e^ ** 3 -^ 1 r m+2ect ,„ +ej3 o r mj „_ e|5 



fc=0 



•5^S , fc_i(0)S fc (S, o )e 8t ' 3 ' f + 1 r m+eai „oT ffl+ea 

a 2 



— £J 2S , ^(9 tfl )r m+ 2e Q ,n+e fl Tm,n~e, 3 + 7i FT7 Tm+e a ,n ° T m+e Q ,ri + 0(B ). 

Cs Qi idi/3,£ + l 



The nullity of the coefficient of £? in this expression, rewritten by using (f76[) , 
leads at once to (I62p . upon doing the relabeling m + e a i— ► to. From it, (|63[) 
follows by duality. Equations — (|6"6")) follow from ([6H)) — (p2l by setting 
/?' ^ /3, a' ^ a and £ = 0. Equations ([57] ) — ([70 ]) follow from (|50] 1 — ([55 ]) by 
setting ft' ^ (3, a' ^ a and forming in each equation the ratio of the cases £ = 
and 1=1, and using the following property of the Hirota symbol, valid for F 
and G depending on a (time-) variable t: 

Equations (fTTj) — (|75|) are just respectively the compatibility equations between 
(fBTj) and JBTJ^^/ , between (1551) and ([B"5jL^w. and between between (|6^|) and 

mi. □ 
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Corollary 6.3 The tau functions T mn and the polynomials P*i^ J (x) — P*i^ n )(x, s* , t*), 
Qm2{y) = Qmn (y, s,t) appearing inY* andY respectively, satisfy the following 
4 formal series identities (6p,p„(n,n*) = (_i)"i+n 2 +-+n <j ,+nl+n5+-+nj„^. 

= 6 p ,0„(n,n*)Tmn(s,t)T m , n ,(s*,t*)(p: { .O(x) Q ( i'J(v) 



V oo 
f3=l 1=0 



m i— > m 

n n — + e^, n* ^ n + ep" — e^ 

S* h-> S a , t/3 i-> */3 - &/3, t* p t-^t l 3 + bp 



+ 



d 2 



Tmn ° Tmn 



= W(«,n*)T ro „(s,t)r m .„.(5*,r)(P i :^', ) ( a; ) QffJd,) 



a=l e=o 



2S'£+25 ctd )r m+e . _ eQ 



772 I > m ea, 777 I ► 777 ~r 6q 
77 I ► 77 e pi , 77* 77 + e^" 

s a i ► s Q a a , s a i ► s a + a a , £^ 
+ 0(a 2 ), 



^"(«»«*)Tmn(s,t)r m . n .(s*,t*) (P*^ m ){x) Q { £2{y) 



m i— > m — e&, m* i— > m + e«j 

77 I ► 77 Cpi, 77* 77 + e/J" 

s* i-> s a , tp^tp- bp, t*ph^tp + bp 



p oo 



X] ^/3^+l^(^) T "»+ea,n+ e( 3 ° 7", 
0=1 £=0 



m — 6a,n- e/3 



+ 0(6 2 



077 d 



/3=1 £=0 



a 2 



t+i 



Tmn ° T mn + 0(b 2 ), 



5p.p„{n,n*)T mn {s,t)Tm*n*{s*,t*) FX'W QS(i/) 



777 I ► 777 

77 77 — Cpi + Bp, 77* 77 + epn — e 

s a i > s Q a Q , s a i » s a + a a , ta i > t/3 



9 oo 



= EE a « 



fi 2 



M-l 



9*3 1 9s a ^+i 



Tmn Tmn + 0(a 2 ), 



a77ii 



" 2 E E a a,e+i s e(dsa ) 

7~m—e a ,n—e0 ° Tm+e Q ,n+e^ + 0(a 2 ). 
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Proof ^From the proof of Theorem 15.11 and (|55|) it follows that 



n i— > n — ep>, n* i— > n* + e^" 
Sq 1 ► s Q Sq, s* i ► Sq. + a c 
tpt->tp — bp, t*p^tp + bp 



= ]T(-l)^ (n) S n} . np+1+k {-2bp)S k {d tf3 )e n( - a ^T m ,^ +ef3 o r„ 

p=l k=0 

and 

9 oo 



^(_l)«r B (m) £ 5 mQ _ m , +1+fc (_2a Q )5 fe (9 s 



„n(o,6) 



ct=l 



fe=0 



and so if we just follow the 4 specializations leading to (|60p - ([55)1 , in order, we 
find the 4 equations of the corollary, in their given order. □ 



7 Examples 

7.1 Biorthogonal polynomials (p = q = 1) 

Given the (not necessarily symmetric) inner product with regard to the weight 
p(x, y) on M. 2 , 

(/0) \g(y)) ■= // f{x)g{y)p{x,y)dxdy 
and the deformed weight 

pt, s (x,y) = e^ (tkyk - Skxk) p(x,y). 

Setting p — q — 1, m = mi, n = ni, with m — n, implies that the indices m, n 
in T mn can be replaced by one single index; namely, set T n : — Tmn-, 

where 



r n (t, s) = det(^x l e~ 



0<i,j<n-l 



Moreover, set = = 1 an( l define the monic polynomials Pn(y) '■= 
Pn\t,s;y) and pn\x) :— pn\t,s;x) (with h~^i the leading coefficient of 
Pn£ l \x)) by 

:= QlnnHy) = y n + --- 

h n -lPn-l( x ) := ^nm H x ) = h n -l x +••• 
The orthogonality conditions Q and (J5|) imply 



pi 1) (y)e Er tkyk ) = for < i < n - 1 



(^Vi 2) (a;)e _l:rsfc:l;fc y j e^° tkyk } = for < j < n - 1 

= 1 for j = n. 
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for all n > 0, from which the bi-orthogonality can be deduced 13 ! 

Pn ] '( x )Pm {y)pt,s(x,y)dxdy = S nm h n . 



From (O, IS]), ® and (JSJ) and from h n — r n+ i/r n , it follows that 

r»(t- [z' 1 ]^) _ (1) 
r n (t,s) ~ Pn [Z) 



K r»(t,s+[z- 1 ]) 

Z -. r 

T n {t, S) 

-n-l T n+l(t + [z _1 ],s) ff p { n\x) 



T n (t,s) JJ M 2 z-y 

-n-l T n+l{t, S - [z- 1 ]) _ ff p { n\y) 



T n (t,s) 

and from ([5T|) , the bilinear identity becomes 



Pt, s {x, y)dxdy 
Pt,s(x,y)dxdy. (77) 



Z — OO 



r„(t, s - [z^ 1 ]) r m (t', s' + [z- 1 ]) e^^~ s >* z m ~ n dz, 

-co 

which characterizes the r- functions for the 2- component KP hierarchy. Equa- 
tions (|77[) and the bilinear identity were obtained in pQ. Indicating the depen- 
dence on t, s in the polynomials, the following inner product can be computed 
in two different ways, leading td 14 l 

T n (t, S )r n+1 (t',s') ([ dxdy pW(t\s';x) P £\t,s;y)eZ?(t*y h -<* k ) p ( x ,y) 



t i— > t — a 
t'i-> t'+a 
s' = s 



2J -2aj + iSj(d t )r n+ 2 o r„ + 0(a 2 ) 



( °° d 2 \ 

= []2^dt^- 1 Tn+lOTn+1+0{a2) )- (78) 

Identifying the coefficients of ctj+i in both expressions and shifting n i— ► n — 1 
yield a first identity; then redoing the calculation above for s i— > s — 6, s' h- > s' + 6 
and f ' = t leads to a second one. All in all we find 

d 2 

Sj(dt)T n+ iOT n -i = -r 2 —— lnr„, 

OSlOtj+l 



13 tt turns out that h n = r n _|_i(t, s)/T n (t, s). 
14 This integral is ^ 0, unless t = t' 
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Sj(d s )T n -x ° r n+ i = ~t^—-— lnr„. 

Specializing the identit y (17311 leads to an identity, which can be expressed as a 
sum of two Wronskianqij and which involves a single tau function: 

f d 2 In t„ d 2 In r» 1 / 9 2 lnr» d 2 In r„ 
\ dtids2 ' dtidsi J t \ dsidt 2 ' dtidsi 

The computation (2.2) was at the origin of the crucial argument (Theorem l5.ip 
in this paper. It illustrates in a simple way what is being done in this paper. 
These equations are used, when computing the PDE for the Dyson, Airy and 
Sine processes ([3]). 




7.2 Orthogonal polynomials 

Given a weight p(z) on K, the symmetric inner product 

if(x) \g(x)) = / f(x)g{x)p(x)dx, 



and the formal deformation by means of an exponential pt(x) :~ p{x)e^i tkZ . 
This is a special case of the previous example, where the deformation only 
depends on i — s; thus t—s can be replaced by t. Then T n (t) is the determinant 
of the moment matrix depending on t = (ti, t<2<, ■ ■ ■), 

T n (t) := detf [ z l+j e^ tkzk p t (z)dz) 

Vr >'• <•.•• " : 

Then, from (|77|). it follows at once that the orthogonal polynomials p n (x) :— 
p n (t;x) are given by 

^nit-lz- 1 ]) 



T n {t) 

- n - 1 T n+1 (t+{z- 1 }) _ f Pn{x) 



= Pn{z) 

p t {x)dx. 



T n (t) 

Moreover, the integral below can be computed in two different ways: on the one 
hand, it is automatically zero, because p n (z) is perpendicular to any polynomial 
of lower degree; on the other hand, for t and t' close to each other, the integral 
can also be developed, using the technique of Proposition 16.21 in t' — t = 2y, 
yielding the following formula 

= T n (t)T n (t') J p n (t; z)p n ^ 1 (t / ,z)p t (z)dz th ^ t _ y 

t' t + y 

15 in terms of the Wronskian {/, g}t = ^j-g — /^f ■ 
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showing that r„ (t) satisfies the KP hierarchy. 

7.3 Orthogonal polynomials on the circle 

Consider the inner product on the circle between analytic functions on S . 

(f(z)\g(z))=<f si ^=-f{z-')g{z) 
and the determinant of moment matrices 



0<k.e<n-l 



T n (t,s) := det^ fc e-^~"« z< 

f . 

g i 27rV-l^ / o<k,e<n-i 



= act f i _J^ z -H£ e En^--«-: 



t r„(i- [z ^s) _ 



Then it follows that 



T n (t,s) 




„n T n(t,S+ [z- 




z — 

T n {t,S) 




n _ 1 T„+i(t+ [Z- 1 ] 




T n (t,s) 




„_ 1 T„+i(i, s- [z~ 




T n (t,s) 





1 27T\/^Tli z - u 



gl 27TIU Z — 

with Pn^(-z) and p„ monic orthogonal polynomials on the circle: 

/ ^-pV(z)p£(z- 1 )=S nm h n , with^ = ^±i. 

The nature of the inner product implies some extra-relationship between the 
orthogonal polynomials 

p^ +l {z)-zp^\z) = p^iO^Hz- 1 ) 
p%{z)-zp^\z) = p%{Q)z n pV{z-'). 
leading to (in the notation of footnote 8) 



{d t )T m +2 O T n (-d s )T m+2 O T r 



T m+2 T n 

In particular, for m = n — 1, 



V h„ J \ h n -ij dti n dsi 
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7.4 Non-intersecting Brownian motions 

Consider N non-intcrsccting Brownian motions X\(t), . . . ,xn(£) in R, leaving 
from distinct points a.\ < ... < and forced to end up at distinct points 
/?!<...< /3jv- From the Karlin-McGregor formula (see [H]), the probability 
that all Xi(t) belong to E C R can be expressed in terms of the Gaussian 
p(t,x,y) = e-^-y^^/V^wi, as follows (0 < t < 1) 



*J (all Xi (t) G B) 
— ^ ' all Xi(t)eE 



(xi(0), . . . ,xjv(0)) = (ai, . . . ,a N ) 
( Xl (l),...,x N {l)) = {(3 u ...,(3 N ) 



1 [ N 

Je n , 



A' 



JJe^-^dxi det 









e * 


det 


e 1 ~ t 




l<i,j<N 





(80) 



l<i,j<N 



The limiting case where several points a and [3 coincide has been the object 
of many interesting studies. It is obtained by taking appropriate limits of the 
formulae above. Just to fix the notation, consider 



(ai,...,a N ) = (ai, ai, . . . , cti, a 2 , a 2 , . . . , a 2 , ■ ■ ■ , a q , a q , . . . , a q ) 

(/3i,...,f3 N ) = (bi,bi,...,bx,b2,b2,...,b2,...,b p ,bp,...,bp), 
where J2l=i a a = Y?p=\ ^3 = and 



«i < 02 < • • • a qi bi < b 2 < . . . < b p , >J m " = n = 

0=1 



0=1 



Then, take the limit of (|8H)l . make a change of variables in the second equality, 
use the standard matrix identity in the third equality 

det K°-« & J>«)i<i,j<„ = det ( a ik)i<i,k<n det ( b ik)i<i,k<n' 

and distribute the integral and the Gaussian over the different columns; this 
yields 



«g (all Xi(t) G E) 

1 /■ *? 

= -TT- \ I ' 2t < 1 - t 'dx z 
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(fa^) , i<mi \ 



x det 



dct 



, fa"*"') 0< 1<m , 

\ 1 < j < N / 

1 r N y 2 

'N JE N M 



< i < 7i! 

1 < 3 < N 



< i < n p , 

1 < i < jv / 



x det 



< i < m t ^ 
1 < j < N 



det 



(//>"-" ) 0< 



/ fa'"") 0<,< ni \ 



• det 



ye 



< i < n p 

1 < j < JV 



\/t(i-*) 
'i-t„ 



N 



dy e 2 y 3 e 



i+3 p (ac+bfi)y 



< i < m c 
< j < tl^ y 1 < a < g 



The numerator of this probability has exactly the form ([T]) evaluated at s c 
tp = 0, with the inner product given by ^ 

(i ! 4W |yV/3(y)} = fdye-^y^e^+'^y, 



upon setting (x) = e a " x , <pp{v) = e?w B.nd dfi(x,y) = S(x— y)e~ y2 I 2 Xe{v) 
By multiplying each of the exponentials e a <* v and e b P y by e~ s< »>* 3 ' arLC j 
e 5Zr '/s.fci/ respectively, it follows that both the numerator and the denominator 
of the probability above, 

T~mn (^1 1 ■ ■ • ) tp ] ^1) • ■ ■ j -^g) 

del i ( / dy e-^-y i+3 e (ha+l ^ v+ ^ (t i ) ' k ~ s "- k)y ) 

I < i < m c 

< j < n„ I 1 < a < q 

1 < <3 < p 



1 < a < q 
1 < /3 < p 



and the same expression for E = R, satisfy the bilinear identity for p + q- 
component KP and, in particular, all the general relations and identities, men- 
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tioned in this paper, namely (fl"2"|) and (jTTj) , (f?2"|) , (f?3")) . Note the equations are 
independent of the set E. 

In particular, for n non-intersecting Brownian motions, departing from the 
origin, with n\ paths forced to end up at —a and paths forced to end up at 
a, we have for < t < 1, 



^(all x(t) £ E) = -^det 



(Mij')o<i<TH— 1, 0<j<ni+7i 2 -l 
{f J >ij)o<i<n 2 -l, 0<j<n 1 +n 2 ~l 



where 

'"ij 

with the change of variables 



fif : = / x i+ ie-^ ±ax dx, (81) 



.2* ~ / 2 
a = a\ and E — EJ— -. 

In a similar way, for several times = to < t\ < . . . < t m < i m +i = 1, 



^"(all Xi{h) EEt,..., all Xi (t m ) e E m ) = ^- det 

An 



I / (My)o<i<ni-l, 0<j<n 1 +"2-l 



(M«)o<»<7»2— 1, 0<j'<rii+n 2 -l 



where 



« m 

/ (x 1 ) j -\x m ) i - 1 e-^T=i^ ±0 ' x rn+T,T=i 1 ^^ + i JJdxt, (82) 



with the change of variables 



2(i m — im-i) „ s / 2(<£+i — ti-i) 

hjc — hip 



(i — t m )(i — t m -i) ' y (t( + i — — tt-i) ' 



and 



, for 1 < j < m — 1. 



y (ij+2 - tj)(tj+i - tj-i) 

We now introduce the inner products 



(/!<?>!= / f(x)g(x)F 1 (x)dx with F!(ar)=e- 

and (m > 2) 



(/l#)m= / /(^i)fl , (^m)-P , m(^i, ■ • ■ ,x m )dxi . . .dx m , 



34 



with 



/ x 2 \ 1 1 

F m (xi, . . . ,x m ) := (TTe - ^ ] e^p^ 1 ^^ 1 c p<! x e x e+i+T,? = 2 E?Li7 r 



The precise form of F m does not matter very much for the purpose of this paper, 
but does play a crucial role in satisfying the Virasoro constraints. 

In these two sets of moments (|8TT) and (|82jl. we insert extra time-parameters, 
as follows, which can then be identified with the moments appearing in ([!]), 



x i+j e - s T±ax±0x 2 e T,T( s k-[l k k 



x 3 e 1 e 



dx 



and 



/ ( \ \ m 



£=1 



x e 



In both cases, we have q = 1 and p = 2, and m := n-i + ri2, leading to the 
introduction of three sets of times Sj := — Su, Mi := — ii$ and Vi := — t2i-Thus, 
from the general theory, the numerator of both probabilities, 



det 



(My)o<i<ni — 1, 0<j<ni+n 2 -l 



(/««) 



0<»<r»2— 1, 0<j<ni+n 2 -l 



satisfies the bilinear identity for the 3-componcnt KP and, in particular, the 
PDE's and thus T raii „ 2 satisfies the single PDE 



9 -. 7Vtl+l,7l2 
■ 111 ■ 



dsi r ni _i !TS2 

A m T "l + 1»"2 

9ui T ni _l )tla 



9 2 l n - 

9si9m 2 mT "i,"2 

d 2 i 
OsiStn Ulr ni,n 2 



• In TV, 



In TV, 



-s-^h — hiT„. „„ dsi tj-S — lnr„. „, 



= 



and the same PDE with Uj replaced by Uj. These PDE's play a crucial role in 
establishing the PDE for the Pearcey process; see [4]. 



35 



X 



The methods developed in this paper should enable one to study more com- 
plicated situations of non-intersecting Brownian motions, as indicated in the 
figure above. The curves in the (x, i)-plane are the boundary of the equilib- 
rium measure as a function of time. When two curves meet, one expects to sec 
a new infinite-dimensional diffusion in that neighborhood, beyond the Pearcey 
process. 
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